Nexus between quantum criticality and the chemical potential 
pinning in high-T c cuprates 
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For strongly correlated electrons the relation between total number of charge carriers n e and 
the chemical potential /i reveals for large Coulomb energy the apparently paradoxical pinning of /1 
within the Mott gap, as observed in high-T c cuprates. By unravelling consequences of the non-trivial 
topology of the charge gauge U(l) group and the associated ground state degeneracy we found a close 
kinship between the pinning of /1 and the zero-temperature divergence of the charge compressibility 
k ~ dn e /dfi, which marks a novel quantum criticality governed by topological charges rather than 
Landau principle of the symmetry breaking. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.72.-h 
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The location of the chemical potential fi is a fundamen- 
tal issue for theories of doped high-T c cuprates. How- 
ever, the understanding of the observed doping depen- 
dence of fj, in cuprates remains a puzzle. In the case 
of LSCO, the angle-resolved photoemission spectroscopy 
studies 0, have shown that in underdoped samples, 
the chemical potential is pinned above the top of the 
lower Hubbard band. The photoemission measurements 
of core levels also shows that /i does not move with hole 
doping in the underdoped region |3(. The observed shift 
in the chemical potential as a function of the electron 
density n e is related to the charge compressibility k or 
the charge susceptibility \c via kT x = (l/n e )(dfi/dn e ) 
and Xc = dn e /d[i, respectively. In Fermi liquid theory 
(FLT) the parameter k is proportional to the effective 
mass m* of the quasiparticles, thus the absence of a shift 
(k — > 00), quite peculiar from the viewpoint of the rigid- 
band model, would imply the divergence of the effective 
mass. The pinning behavior of (i can be made plausible 
for a system where charge carriers are segregated, e.g., in 
a stripe form However, it is hard to imagine how the 
property m* — > 00 might arise due to inhomogeneities, 
since the mass enhancement appears to result from the 
physics of the uniform liquid 5]. Anyway, the behavior 
of y, is quite peculiar from the viewpoint of the FLT of 
the metallic states and signals a dramatic reorganization 
of the electronic structure of cuprates with doping. 

Given the fact that the electron density appears (be- 
sides the temperature T) as a control parameter, it fol- 
lows from scaling and thermodynamical considerations 
that the existence of quantum phase transitions and 
associated zero-temperature critical points might be a 
part of the physics of high-T c cuprates. Indeed, in 
cuprates there is clear evidence for the existence of a spe- 
cial doping point in the lightly-overdoped region where 
superconductivity is most robust. Such behaviour indi- 
cates this point could be a quantum critical point (QCP) 
while the associated critical fluctuations might be respon- 
sible for the unconventional normal state behaviour 



However, it is unclear whether this QCP is "truly crit- 
ical" in the sense that it is characterized by universal- 
ity and hyperscaling. The resemblance to a conventional 
QCP is hampered by the lack of any clear signature of 
thermodynamic critical behavior. Experiments appear to 
exclude any broken symmetry around this point although 
a sharp change in transport properties is observed ||- 

Understanding these novel phenomena in cuprates are 
challenging tasks. In the celebrated Hubbard model, 
which epitomizes the character of strongly correlated 
electrons the degeneracy of energy states renders the 
usual perturbation theory inapplicable. However, there is 
an intimate relation between the ground state degeneracy 
and the topology of the system- much of the quantum in- 
formation is encoded in phase relations of the many-body 
wave function The non-integrable Berry phase fac- 
tors of the topological origin that the eigenvectors of 
the Hubbard Hamiltonian possesses can produce the in- 
terference akin to the Aharonov-Bohm effect [hJ in an 
electromagnetic field giving rise to more exotic elemen- 
tary excitations 

In this work, we explore Mott transitions from the non- 
magnetic insulator to a superconductor induced by dop- 
ing and show that the process is governed by the topo- 
logical structure of the electromagnetic compact gauge 
U(l) group. As a result collective instanton excitations 
of the phase field (dual to the charge) arise with a great 
degree of stability, governed by gauge flux changes by 
an integer multiples of 2ir, which labels the ground state 
degeneracy. The associated abrupt transition between 
differnt "vacua" allows us to make link between the un- 
usual behavior of the chemical potential and a novel type 
of quantum criticality that goes beyond the paradigm of 
the symmetry breaking. 

We consider an effective one-band electronic Hamil- 
tonian on a tetragonal lattice that emphasises strong 
anisotropy and the presence of a layered CuC"2 stacking 
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sequence in cuprates: H — Ht-j + TLu + H±, where 
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Here (r,r'} and ((r, r')) denotes summation over the 
nearest-neighbour and next-nearest-neighbour sites la- 
belled by 1 < r < N within the CuO plane, re- 
spectively, with t, t' being the bare hopping integrals 
t' > 0, while 1 < £ < N± labels copper-oxide lay- 
ers. The operator c^(r)(c Q ^(r)) creates (annihilates) 
an electron of spin a at the lattice site (r, £), S$(r) = 

S Q /3 c L( r ) cr Q/3 c ^( r ) ( a = x ->ViZ) stands for spin and 
n^(r) = nie(r) + n^r) number operators, respec- 
tively, where rw(r) = c^(r)c Q ^(r), [i is the chem- 
ical potential and J the antifferomagnetic (AF) ex- 
change. The Hubbard term is Jiu = J2 er Un^e(r)nie(r) 
with the on-site repulsion Coulomb energy U, while 
H± = — J2rr' ^J-( rr ') c i^( r ) c af+i( r ') facilitates the inter- 
layer coupling, where t± is the interlayer hopping with 
the c— axis dispersion e±(k,k z ) = 2t±(k) cos(ck z ), while 
ij_(k) = t± [cos(ak x ) — cos(ak y )) 2 as predicted on the ba- 
sis of band calculations 0] . We write the partition func- 
tion Z = J [DcDc] e - 5[5 ' c] with the action 
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Jo „ 



c a t{vT)d T c a i{vT) +H(t)} (2) 



using coherent-state fermionic path integral over Grass- 
mann fields depending on the "imaginary time" < 
r < j3 = 1/fcsT. Further, we write the Hub- 
bard term in a SU(2) invariant way as T~Lu{t) = 
^E rf {(l/4)«£ 2 (rr) - [^(rr).S,(rr)] 2 } with charge- 
U(l) and spin-SU(2)/U(l) sectors, where the unit vector 
£lf(rr) sets varying in space-time spin quantization axis 
yj| . In the following we fix our attention on the U(l) 
invariant charge sector [i ll . We consider now the follow- 
ing resolution of unity using the Fadeev-Popov method 



0, 1 ee f[DQ]S[Q-n] = f[W] e iX*IfW-»)v, 
where n(rr) = ^2 a c Q (rr)c a (rT) represents the particle 
number, Q(yt) is the collective variable and iV(rr) is 
fluctuating (in space and time) imaginary chemical po- 
tential conjugate to the local particle number niivr) . 
The field Vi(tt) can be written as a sum of a static 
Voi(r) and periodic function V(rr) — Vo(r) + V(rr) us- 
ing Fourier series V(rr) = {1/ P)Y^=i\V{ YUJ n)e l ^ n ' r + 
c.c] with u> n = 2Trn/(3 (n — 0, ±1,±2) being the 
(Bose) Matsubara frequencies. Now, we introduce the 
phase (or "flux") field <^(rr) via the Faraday-type re- 
lation ^(rr) = 8i M rT ) = Vf(rr) to remove the imag- 



inary term i J Q dr^(rr)n^ (rr) = i J Q drVg (rr)n^(rr) 
for all the Fourier modes of the V^(rr) field, except 
for lo„ = by performing the gauge transformation to 
the new fermionic variables f a i(rr), where c a ^(rr) — 
e «<Mrr) J Q ^(rr), which indicates that the electron ac quir e 
a phase shift similar to that in the electric AB effect . 
The electromagnetic U(l) group governing the phase field 
is compact, i.e. 4>i{yt) has the topology of a circle (Si), 
so that topological effects can arise due to non-homotopic 
mappings of the configuration space onto the gauge group 
S\ — ► U(l). To explicate the composite nature of the 
physical electron field we write explicitly: 



c a e{ rT ) = ex P 



I / dT'V Dl (vT') 




e <7« ("•)/„,(„-). (3) 



The first term in the exponential in Eq.Q is the usual 
dynamical phase factor where 6g(rr) = Vd(l{yt) and 
0i(y(3) — 9i(r0) . The second one is the non-integrable 
Berry phase factor |9(: 7fs(rT) = 2irTme(r)/ f3, where 
mi(r) marks the integer U(l) winding number. Work- 
ing backwards we can now formulate the the path inte- 
gral for the partition function. For this we concentrate 
on closed paths (or world lines)in the imaginary time r 
that start at position r£ at imaginary time r = and 
end at the same position at r = j3, which fall into dis- 
tinct, disconnected (homotopy) classes labelled by the 
winding number |Tjj. Homotopically distinct paths can 
be summed according to various possibilities for inequiv- 
alent quantizations (superselection sectors) according to 
the formula: 
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The first order in time derivative term in Eq.Q) is just 
the topological action which does not affect the equa- 
tion of motions but influences statistics [l2(. The gauge 
transformation in Eq.© introduces phase factor into 
the hopping elements of the hamiltonian which frustrate 
the motion of the fermionic subsystem. However, when 
charge fluctuations become phase coherent is signalled by 
(e J< ^( rT )) ^= the frustration of the kinetic energy is re- 
leased. To proceed, we trace over the fermionic degrees 
of freedom in Eq. J2J and introduce the unimodular com- 
plex scalar ze(rr) — e i ^^ rT - ) via another Fadeev-Popov 
resolution of unity 
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where the unimodularity constraint can be 
imposed with a real Lagrange multiplier A. 
The partition function then becomes Z = 
J[V 2 z]U re S(\z e (rT)\ 2 -l)e- s ^ z ' z *\ where S[z,z*] = 

£ qw „ z q( w '») r q 1 ( w «)^( w '')' Tlie P hase " 

' Eq^„ r qK) 



PNNI ^qu„ z q( w n)r q (w„)^ q (o; n ). 

coherence boundary is 1 = 



PNN± 



Here, q = (k, k z ) and at criticality, the condition 
r q =o(k->n = 0)|a=a c = fixes the Lagrange parameter 
while 

TqVn) = A-S(q,w„)+ 7 - 1 K) 
E(q, uj n ) = J|[ (A) cos(afc x ) cos(ak y ) 

+ E^^[4(A) + Ji(A)™(*)] ,(6) 
q' 

where £k = cos(ak x ) + cos(ak y ) Furthermore, 7o(w„) 
is the Fourier transform of the bare phase correlator 
( e -*[0i(n-)-0 4 /(rV)]^ Q originating from the kinetic and 
topological part of the action in Eq. (0} 
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where Z = exp(-^-)0 3 (|^, e^/ 4 ) and 3 ( U)9 ) = 
g m2 e 27rm ™ is the Jacobi theta function [3 



E + oo 
m— — oo ' 

which is /3-periodic in the imaginary-time r as well as in 
the variable 2fx/U with the period of unity. The micro- 
scopic phase stiffnesses to the lowest order in the hopping 
amplitudes are given by 
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FIG. 1: (Color online) The critical temperature T c as a func- 
tion of the chemical potential /i together with the density 
plot of the charge susceptibility \ c = Uxc/2 for t* = 0.5eV, 
t'*/t* = 0.3, t* ± = O.OleV, J = 0.15eV and V = 4eV. Inset: 
Xc and the occupation number rib for T = 0.1C7. 



where p = fj. — rifU/2, rif = (/ Q (rr)/ Q ^(rr)) is the oc- 
cupation number for /— fermions and v n = n(2n + 1)//? 
(n = 0, ±1, ±2) stand for the Fermi Matsubara frequency. 
The stiffnesses in Eq. (jHJ) rest on the "<i-wave" pair ampli- 
tude A(k) = |A|?7k where rj^ — cos(ak x ) —cos(ak y ) due to 
the in-plane momentum space pairing of the /-fermions. 
A Gorkov-type decoupling of the AF exchange term in 
Eq. ||TJ using the valence bond operator readily gives 

for the gap parameter: 1 = ^ k tanh []|/3£k] with 
the quasiparticle spectrum, Sy = [e|(k) — p] 2 + |A(k)| 2 . 
Here, e|(k) is effective in-plane band dispersion narrowed 
due to the frustrated motion of the carriers in the fluc- 
tuating "bath" of U(l) gauge potentials, so that the ac- 
tual tight-binding parameters are "dressed" ones t* x = 
tx ^[Mrr)~Mr'r)]^ where ty . = t,t',t x are the bare 
band parameters. Now, by differentiating the free energy 
before after implementing of the gauge transformation 
in Eq.(|3J) we obtain for the electron occupation number 
n e = (c ct (rr)c ct ^(rr)) the result n e = n/ + 7i(, — 2fi/U. 
Here, rib — 2[i/U+ (2/iU){4>{tt)) is the occupation num- 
ber given by the mean value of the winding numbers 
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FIG. 2: (Color online) The chemical potential fi as a function 
of the occupation number rif, for various temperatures T as 
indicated in the plot and model's parameters as in Fig.l. The 
value of p stays within the charge gap as ni changes. Inset: 
the density plot of rib as a function of fi and T. 



ior of /i characteristic for a bosonic system. 

To conclude, new type of quantum numbers must be 
invoked to explain topologically induced QCP in cuprates 
and the associated pinning of the chemical potential. 
Topological effects arise as stable, non-perturbative, col- 
lective excitations of the phase field (dual to the charge), 
which carry novel topological characteristics. These 
are the winding numbers of U(l) group: w^(r) = 
la drifie (rr) that become topological conserved quan- 
tities. It is exactly the appearance of these topologi- 
cal charges that render the system "protected" against 
small changes of the hamiltonian's parameters. This 
novel conservation does not arise just out of a symme- 
try of the theory (as "conventional" conservation laws 
based on Noether's theorem) but it is a consequence of 
the connectedness, i.e. topology of the phase space, re- 
lated to the topological properties of the associated sym- 
metry group. 

This work was supported by the Polish Science Com- 
mittee (KBN) under the grant No. 2P03B 009 25 in years 
2003-2005. 



where we made use of the relation 
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In the limit of strong (weak) correlations n e interpo- 
lates between topological rib (fermionic nf) occupation 
numbers. In the large-?/ limit ix — ► ntU/2, so that 
n e — * nb and the system behaves as governed entirely by 
U(l) topological charges which play the role of "quasi- 
particles" . It is straightforward now to calculate the 
charge compressibility k. The result is given in Fig.l 
along with the outcome for the superconducting phase 
boundary. We see the evolution of K with decreasing 
n e , (i.e. hole doping) from the Mott insulator |2jj with 
k = (at 2fj,/U = 1) to a point of degeneracy on the 
brink of the particle occupation change at 2fi/U = .5 
where k — oo at T = 0. This is also the point on the 
phase diagram from which the superconducting lobe em- 
anates. It is clear that, the nature of the divergence of 
k here has nothing to do with singular fluctuations due 
to spontaneous symmetry breaking as in the "conven- 
tional" phase transition. Rather, this divergent response 
appears as a kind of topological protection built in the 
system against the small changes of /i. Further, k — > oo 
implies that the and dii/dn e becomes vanishingly small 
at T = which results in the chemical potential pinning, 
see Fig. 2, where the temperature dependence of \i is also 
shown in the inset. In cuprates fi(T) can b accessed via 
the measurement of the work function [2lJ : the available 
data for YBCO clearly indicates the temperature behav- 
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